Frustration-induced non-Fermi-liquid behavior in a three-impurity Kendo model 
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A Kondo model for three spin-i impurities placed at the vertices of an equilateral triangle is studied 
using the numerical renormalization-group method. The impurity spins can form two frustrated 
doublets which, because of the rotation symmetry, couple equally to the conduction band. Over a 
broad range of antiferromagnetic RKKY interaction strengths, this degeneracy drives the system to 
one of two novel, non- Fermi- liquid fixed points. Potential scattering is a relevant perturbation at 
one of these fixed points, but is marginal throughout the second regime. 



> 
O 

O 

\o 



o 
o 



75.20.Hr 

Much of the rich behavior of heavy- fermion compounds 
is thought to stem from competition between order- 
ing of local moments via the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interaction, and quenching of those mo- 
ments due to the Kondo effect. Theoretical understand- 
ing of the subject is far from complete, ft is widely be- 
lieved that the essential physics is contained in the pe- 
riodic Anderson Hamiltonian, but this model has been 
solved only at the mean- field level Insight has also 
been gained from exact solutions for one or two mag- 
netic impurities in simple metals |^], and for lattices in 
the limit of infinite spatial dimensionality Q]. 

Recently, attention has focused on rare-earth and ac- 
tinide systems which exhibit a specific heat coefficient 
C/T that diverges weakly with decreasing temperature 
T — rather than remaining constant, as in a (heavy) 
Fermi liquid — and a resistivity that is linear — in- 
stead of quadratic — in T Explanations proposed 
for these non- Fermi- liquid (NFL) properties include: (1) 
the single-impurity Kondo effect with a disorder-induced 
distribution of energy scales (2) isolated magnetic im- 
purities coupled to two or more conduction bands [Q ; and 
(3) proximity to a T = magnetic-ordering transition Q . 

A fourth possibility, overlapping with but not identi- 
cal to (3) above, is that lattice effects, such as period- 
icity or geometric frustration, may induce NFL behav- 
ior where none exists in the one-impurity limit. Few- 
impurity systems may provide clues as to the plausibility 
of this scenario. The two-impurity Kondo problem yields 
NFL physics — but only at particle-hole symmetry and 
with fine-tuning of the RKKY coupling [^,0 . 

In this paper, a more robust form of NFL behavior 
is identified in a model for a single conduction band in- 
teracting with three identical spin-i impurities arranged 
in an equilateral triangle. Magnetic frustration prevents 
even strong antiferromagnetic RKKY interactions from 
locking the impurities into a total spin singlet; instead, 
four degenerate spin-i configurations remain. Using the 
numerical renormalization-group method, we find that 
over a wide range of RKKY couplings, Kondo screening 
of these impurity states is governed by two NFL fixed 
points. Particle- hole asymmetry destabilizes one of the 



regimes, but is a marginal perturbation in the other. Our 
results suggest that spatial symmetry might give rise to 
NFL behavior in more realistic models without the need 
for parameter tuning. (Similar considerations have led 
to the study of electron-assisted tunneling between three 
atomic sites [0 as a possible explanation for the anoma- 
lous conductance of certain point contacts fl^.) 

We start with a Hamiltonian, H = Hband + Hint, de- 
scribing a non-interacting conduction band coupled lo- 
cally to three spin-i impurities. Si {i = 0,1,2), which 
lie at the vertices of an equilateral triangle. The band 
is assumed to be isotropic in momentum space, and to 
extend in energy over a range ±D about the Fermi level 
(taken to be energy e = 0). The interaction term is 
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[Vn{vi) - Js(r,) • S,] 
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where n{ri) [s(ri)] is the electron number [spin] at impu- 
rity site i. The exchange coupling J is assumed negative, 
whereas the potential scattering V can be of either sign. 

Symmetry permits the impurities to couple to just six 
orthogonal, spatially localized conduction states, annihi- 
lated by operators Cha, where cr = ±i labels the spin z 
component, and h £ {0, 1, 2} specifies the "helicity." He- 
licity is analogous to parity in a two-impurity problem, 
i.e., a helicity-/i wave function acquires a multiplicative 
factor e*^'^'*/'^ under a rotation of 27r/3 about the cen- 
ter of symmetry. The combined states of three spin-i 
impurities decompose into a quartet having total spin | 
and helicity h — Q, plus two doublets of total spin i. 
The doublets can be constructed so that one has h = 1 
and the other h = 2, in which case the Hamiltonian con- 
serves total helicity (modulo 3), and is invariant under 
interchange of helicity labels I and 2. 

Applying these observations, Eq. (j^) can be rewritten 

H^nt = Vnun + Vi{ni + na) - Jqo Sqo-So 

- Jll(sii -I- S22)-So - Jl2(si2-Si -I- 821-82) 
--'"oi[(soi + S20)-Sl + (Sio + S02)-S2], (2) 



where rih = J2o 



is the number of helicity-/i elec- 



trons at the impurity sites, Shh' = J2ctct' ^hah^ 
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gives the spin of an electron scattering from helicity h' 
to h, and S^h = X]j=o e'^'^-^'^'^Sj is the spin of the im- 
purities as their heUcity is raised by A/i. The couphngs 



Vh = VFl, Jhh' = JFhFh 



(3) 



depend on the band dispersion A:(e), the density of states 
p(e), and the impurity separation i?, through form factors 
K = \S depie) [1 + (35^,0 - 1) sin(A:i?)/(fci?)] . 

This problem can be solved using a generalization of 
the non-perturbative renormalization-group method de- 



veloped for the one- impurity Kondo model 1 13|. The con- 
tinuum of electronic energies is replaced by a discrete 
set, ibDA-", n = 0,1,..., where A > 1. The Hamil- 
tonian is then mapped, via the Lanczos method, onto 
a tight-binding model for a one-dimensional chain with 
the impurities coupled to one end. The discretization in- 
troduces a separation of energy scales which allows the 
Hamiltonian to be solved iteratively, starting with a sin- 
gle site. At each subsequent step, one site is added to the 
chain, while the effective temperature T falls by a factor 
of A^/^. This continues until the energy spectrum be- 
comes scale-invariant, i.e., a stable fixed point is reached. 

The discretized Hamiltonian contains an infinite set of 
parameters which depend on the band shape and the im- 
purity separation R. We argue, along lines advanced for 
the two-impurity Kondo problem [l4| , that these band 
parameters can be replaced by standard values calcu- 
lated for p{e) = constant and R — oo. All possible 
low-temperature regimes of the original model can be 
explored provided that one treats the couplings entering 
Eq. (|^) as independent parameters, freed from the con- 
straints implied by Eqs. (^). In particular, Jn and J12, 
although equal as originally defined, should be regarded 
as separate couplings because no symmetry of the prob- 
lem prevents them from renormalizing independently. 

We have solved the discretized version of this problem 
for different initial values of the couplings Vb, Vi, Joo, 
Jii, J12, and Joi. Even taking advantage of conserved 
quantities to break the Hamiltonian matrix into block- 
diagonal form and thereby to reduce CPU time, we have 
generally found it practical to retain no more than 1200 
states after each iteration. This corresponds to keep- 
ing all energies E < 5T, which allows identification of 
fixed points, but is insufficient for direct computation of 
physical properties as thermodynamic traces. The good 
quantum numbers are total helicity, H; total spin and its 
z component, S and Sz', and total "isospin" ||lO| or "ax- 
ial charge" Q and its z component, / and Iz ]15|]. Total 
isospin is only conserved at particle-hole symmetry, so 
most runs were performed with Vo = Vi = in order to 
benefit from the additional quantum number. 

The three-impurity Kondo problem has a rich phase 
diagram, as might be expected given the large num- 
ber of independent parameters. In this paper, we focus 
on the stable low-temperature regimes of the particle- 
hole-symmetric model. Figure ^ shows a schematic 
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FIG. 1. Schematic phase diagram for the three-impurity 
Kondo model at particle-hole symmetry. Here if is a mea- 
sure of the RKKY interaction, and a measures the relative 
strengths of the couplings Joi and J12 in Eq. (^). Circles rep- 
resent stable fixed points, and dashed lines mark the borders 
between difi'erent regimes of low-temperature behavior. 



phase diagram for this case, projected onto the plane 
spanned by two dimensionless combinations of the four 
exchange couplings entering Eq. (|^): The first, K — 

("^00 + ^t^ii — J'\2 ~ 2>^oi)/(>^oo + 2J^]^ -l- 2J^2 + 4Jq]^), 
measures the strength of the RKKY interaction. As 
K ranges from — ^ through to +1, the static impu- 
rity correlation fimction (i Xli<j ' ^i) varies smoothly 
from —1/4 through to -1-1/4. The second combina- 
tion, a = (Joi — ^i2)/(>/oi + >/i2), measures the difference 
between the two couplings that dominate the region of 
antiferromagnetic RKKY interactions, K 

There are three distinct regions in Fig. |l|: 

Independent- Kondo regime: For ferromagnetic and 
weakly antiferromagnetic RKKY couplings, the system 
renormalizes to a stable fixed point at which interac- 
tions between the impurities are essentially irrelevant. 
The ground state is a spin singlet, and all electrons at 
the Fermi level undergo a phase shift Shi^ = 0) = 7r/2. 
Thus, the low-energy excitations are those for three im- 
purities at infinite separation, each independently under- 
going a Kondo effect. Particle-hole asymmetry is exactly 
marginal in this regime, just as it is in the one-impurity 
Kondo model. 

Frustrated-Kondo regime: In a region of medium-to- 
strong antiferromagnetic RKKY interactions where Jqi is 
the dominant coupling, the system flows to a novel, stable 
NFL fixed point at which electrons of all three helicities 
screen the frustrated spin-i impurity configurations. 

Table Q gives the finite-size spectrum at this fixed point, 
computed for A = 3, a value that represents a near- 
optimal compromise between discretization errors (which 
grow with A) and numerical rounding errors (which be- 
come larger as A approaches unity) |l^. The energies 
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TABLE I. Low-lying states from the finite-size spectrum 
at the frustrated-Kondo fixed point. Energies E axe given for 
odd- and even-length chains. States are labeled by their spin 
S, isospin /, and helicity H\ an asterisk stands for one state 
each of helicity 1 and 2. See the text for further details. 



have been multiplied by a factor of 0.625, chosen to map 
the smallest splitting in the spectrum of the free conduc- 
tion band (without any impurity) onto its value in the 
continuum limit, A ^ 1. The scaled energies E are ex- 
pressed in units of Trvj^/L, where L is the number of sites 
in the chain, and vp is the Fermi velocity. 

The regularity of the levels in Table | is reminiscent 
of a Fermi liquid, but several details — particularly the 
absence of states in the odd-L spectrum at i? « 0.4, and 
the ordering of the even-L states for E < 0.7 — prove 
impossible to explain within this framework. 

The rate at which the eigenenergies approach their 
fixed-point values also points to an NFL description: By 
examining the L-dependence of the spectrum, and using 
the relation T DA^^/'^ |l|], we deduce that the lead- 
ing irrelevant operator at this fixed point varies like T^, 
with A = 0.21 ±0.03 (c.f. A = 1 for a Fermi liquid). This 
exponent should be reflected in anomalous temperature 
dependences of physical properties. 

An NFL interpretation of the frustrated-Kondo regime 
is supported by evidence that the fixed point occurs at 
finite, non-zero exchange couplings. Empirically, one 
can eliminate all spin-| impurity configurations, and set 
Joo = Jii = J12 = in Eq. (^), to arrive at a one- 
parameter model which, for all —D < Joi < +D (at 
least), reproduces the fixed-point spectrum in Table ^ It 
can be shown using poor-man's scaling fl^ that the Jqq 
and Jii terms in Eq. (^ are regenerated under renormal- 
ization; however, J12 remains identically zero. In princi- 
ple, one should be able to locate the fixed point by ad- 
justing Joo, Jii, and Joi until the energy spectrum does 



not change at all as the tight-binding chain is extended. 
In practice, the band discretization introduces irrelevant 
operators which renormalize the spectrum at small val- 
ues of L. We therefore associate the fixed point with 
the couplings which yield the smallest deviation from the 
asymptotic spectrum at some standard chain length, usu- 
ally L — 20. This procedure gives Jqq = ^01 = — 0.8D 
and J*i = +OAD (all to within ±0.1D), and J*2 = 0. 

Throughout the frustrated-Kondo regime, potential 
scattering is found to act as a marginal perturbation 
about the particle-hole symmetric limit, shifting the en- 
ergy of each state by an amount proportional to its charge 
quantum number, J^. Thus, for moderate values of Vq 
and Vi [see Eq. (||)], this regime retains its NFL char- 
acter. Very strong potential scattering suppresses the 
Kondo interaction, however, and drives the system to a 
weak coupling Fermi liquid. The crossover between these 
extremes is hard to analyze given the limited number of 
states we can retain in our calculations, but this issue 
merits future attention. 

The evidence presented above demonstrates that the 
frustrated-Kondo regime has an NFL fixed point which 
is stable with respect to the RKKY interaction and 
is marginal under particle-hole asymmetry. (However, 
breaking the real-space symmetry of the impurity loca- 
tions will lead to the recovery of Fermi-liquid behavior be- 
low some crossover temperature.) This fixed point does 
not seem to map onto any fixed point known from pre- 
vious studies of impurity models. One likely candidate, 
the overscreened limit of the 3-channel Kondo model for a 
spin-i impurity, is ruled out because it predicts a ground 
state spin S = 1 for L odd, and S = for L even. 

Isospin-Kondo regime: In a second region of the 
parameter space — characterized by antiferromagnetic 
RKKY interactions, but dominated by J12 rather than 
Joi — the system flows to a stable NFL fixed point which 
is related, by interchange of spin and isospin degrees of 
freedom, to the one-impurity, two-channel Kondo model. 

Table |ll| gives the finite-size spectrum at this fixed 
point, once again computed for A = 3 and scaled by 
a factor of 0.625. The energies and quantum numbers 
of all the low-lying states can be reproduced by convo- 
luting the free-fermion spectrum for decoupled helicity-0 
electrons with a spectrum derived from the NFL fixed 
point of the one-impurity, two-channel model [p^ : For 
each state at the two-channel fixed point, labeled by spin 
J, flavor Jf and charge Q [p7| , p^ , one constructs a three- 
impurity state with the same energy, having quantum 
numbers S = Jj: I = S; and iJ = 0, 1 or 2, depending 
on whether Q = 0, Q > or Q < 0. 

Insight into these observations can be gained from an- 
other one-parameter model. Spin-| impurity configu- 
rations are eliminated, just as for the frustrated-Kondo 
fixed point, but here J12 is the only non-zero coupling. 
For all —D < J12 < +D, this model flows to the fixed- 
point listed in Table ||. Since Joo = Joi = 0, the helicity- 
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stronger antiferromagnetic RKKY couplings, the three- 
fold spatial symmetry of the model leads to two novel 
low-temperature regimes of non-Fermi-liquid behavior, 
one of which is marginally stable against particle-hole 
asymmetry. Further work is required to identify the phys- 
ical properties of the stable non-Fermi-liquid regime, and 
to determine whether or not it can be reached starting 
from any realistic band structure. 
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TABLE II. Low-lying states at the isospin-Kondo fixed 
point. See Table ^ for an explanation of the entries. 

electrons must be decoupled at the fixed point. The in- 
teraction between the remaining degrees of freedom can 
be written 

1 ^ ^ 

/^,i/=l h,h' = l a, a' 

where S and H are commuting SU(2) operators acting on 
the impurity spin and helicity, respectively. An isotropic 
{Jul/ = J) version of Eq. (^ has been considered in com- 
bination with pure spin and flavor interactions in a gen- 
eralized two-channel model [|l^ , but we are aware of no 
study of the case J^i, = J(l — Si,,^) which arises here. 

Unlike the other two fixed points, the isospin-Kondo 
fixed point is unstable with respect to particle-hole asym- 
metry. Specifically, Vi entering Eq. (||) plays the role of 
the magnetic field in the standard two-channel Kondo 
model, driving the system to a Fermi-liquid state. Since 
there is no reason to expect Vi to be zero in any real 
system, this instability implies that the isospin-Kondo 
regime is unlikely to be observed in practice. 

In the preceding discussion, we have not specified the 
precise positions of the boundaries between the three 
regimes shown in Fig. |^. These depend in a fairly compli- 
cated fashion on the two combinations of exchange cou- 
plings orthogonal to K and a. We note, though, that the 
frustrated-Kondo and isospin-Kondo regimes expand as 
the Kondo coupling J decreases in magnitude, i.e., as the 
RKKY coupling grows to dominate the single-impurity 
Kondo effect. We suspect that in the limit J ^ 0, each 
NFL regime occupies an entire quadrant of Fig. |^. 

We have presented numerical renormalization-group 
results for a model of three magnetic impurities interact- 
ing with a non-magnetic metal. For ferromagnetic and 
weakly antiferromagnetic RKKY interactions, the impu- 
rities undergo essentially independent Kondo effects. For 
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